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Abstract 



We demonstrate an unambiguous and robust method for computing fermionic 
corrections to the energies of classical background field configurations. We 
consider the particular case of a sequence of background field configurations 
that interpolates continuously between the trivial vacuum and a widely sep- 
arated soliton/antisoliton pair in 1+1 dimensions. Working in the contin- 
uum, we use phase shifts, the Born approximation, and Levinson's theorem 
to avoid ambiguities of renormalization procedure and boundary conditions. 
We carry out the calculation analytically at both ends of the interpolation 
and numerically in between, and show how the relevant physical quantities 
vary continuously. In the process, we elucidate properties of the fermionic 
phase shifts and zero modes. 
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I. INTRODUCTION 



1+1 dimensional models have pointed out many subtleties in the computation of one-loop 
fermionic quantum corrections to the energies of classical field configurations. Many works, 
usually in the context of supersymmetric theories [1], have reached conflicting conclusions 
about the correct boundary conditions, regularization and renormalization procedures, and 
bound state and continuum contributions. Since the one-loop energy is given by the cancel- 
lation of a divergent sum of zero-point energies against a divergent counterterm, one must 
be extremely careful to fix the counterterm with a definite renormalization scheme in order 
to obtain a physically meaningful result. 

In this paper we extend the methods of Ref. [2,3] to this case. We again find that 
1 + 1 dimensional models provide an excellent testing ground for our methods, and can give 
considerable insight into the causes and resolutions of the disagreements among the previous 
results. As in our previous works, we will use a continuum formalism based on phase shifts, 
their Born approximations, and Levinson's theorem. Relying on such physical quantities 
results in a procedure that is both analytically unambiguous and numerically robust. As a 
result, we will be able to compare numerical results in generic, but not exactly soluble cases 
to analytic results in specific, exactly soluble cases. 

We will consider a Majorana fermion coupled to a real scalar background field, with 
a definite choice of bosonic and fermionic potentials, though our method applies equally 
well for any other choice of potentials, as well as for complex scalars and Dirac fermions 
(including cases without C or CP invariance). We will consider a continuous sequence 
of scalar field configurations that interpolates between a trivial background and a widely 
separated soliton/antisoliton pair, always with the same, trivial field configuration as x — > 
±oo. In this way we can unambiguously track physical quantities such as phase shifts and 
bound state energies, staying in the trivial sector of the theory throughout the process. 
In so doing, we obtain results numerically that continuously approach the results obtained 
analytically in Ref. [4]. 

In §11 we develop the methods we will need to compute phase shifts, which in turn allow 
us to compute the renormalized one-loop energy as a sum over zero-point energies —\oj- In 
§111 we apply these methods to the specific sequence of field configurations we have chosen 
and obtain numerical and analytical results. In §IV we summarize our results, and connect 
them to previous and future works. The Appendix contains techniques parallel to those of 
§11 but for a single soliton instead of a soliton/antisoliton pair. 

II. FERMIONIC SMALL OSCILLATIONS AND PHASE SHIFTS 

We consider a Majorana fermion \& interacting with a scalar background field 0, with 
the classical Lagrangian density 



where is the Lagrangian density for the <fi background field. As discussed in the Intro- 
duction, we can treat any bosonic potential. To be specific, we choose 
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which has classical soliton (antisoliton) solutions satisfying 



(f)'(x) 
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giving 



(f){x) 
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For the most part, we will be treating as a classical background, so that its dynamics 
will be unimportant. However, we note that choosing the Lagrangian density of Eq. (2) 
causes the bosonic and fermionic degrees of freedom to be related by supersymmetry for a 
background field that is a solution to the equations of motion (such as a single soliton or an 
infinitely separated soliton/antisoliton pair). In what follows, we will find it instructive to 
compare the bosonic and fermionic small oscillations spectra. The supersymmetric model is 
discussed in more detail in Ref. [4]. 

The one-loop corrections to the energy due to fermionic fluctuations will be given by 
the appropriately renormalized sum of the zero-point energies, — of the fermionic small 
oscillations. The spectrum of fermionic small oscillations in a background 4>o(x) is given by 
the Dirac equation, 



labeling the scattering states. 

We will choose the convention 7 = a 2 and 7 1 = W3 so that the Majorana condition 
becomes simply = We note that since our Lagrangian is CP invariant, all of our 
results for the spectrum of a Majorana fermion can be extended to a Dirac fermion simply 
by doubling. 

Our primary tool for analyzing the small oscillations will be the phase shifts. We can 
solve for the phase shifts of any fermionic potential Vf(x) that satisfies Vf(x) = Vp(— x) 
and Vp{x) — > m as x — > ±00 by generalizing the methods of Ref. [2] to fermions. This form 
will be useful for considering our example of a sequence of background field configurations 
that continuously interpolates between the trivial vacuum and a widely separated soliton- 
antisoliton pair. In the limit of infinite separation, the phase shift for the pair goes to twice 
the result for a single soliton. For comparison, we also do the computation for a single 
soliton directly in the Appendix. 

We define the parity operator acting on fermionic states as P = II7 , where II sends 
x — > —x. P commutes with the Hamiltonian, so parity is a good quantum number. As a 
result, we can separate the small oscillations into positive and negative parity channels, now 
restricted to k > in the continuum. 

We parameterize the fermion wavefunction by 




(5) 



where 




V 'uS 1 — m 2 is the momentum 




(6) 
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where 9 = tan 1 ^ and v and ( are complex functions of x. We then find the phase shift in 
each channel 5 ± (k) by solving Eq. (5) subject to the boundary conditions 



We obtain 

^(fc) = -Re „(0) + °- + 1 log 

r(fc) = -itei/(0) + ? + itogi^ (8) 

where F = ^(V F {0) - ik + ii/(0)*) and satisfies 

+ v\xf + 2fez/(a;) + V F {x) 2 - Vf(x)' - m 2 = (9) 

with the boundary condition that v{x) and v'{x) vanish at infinity. The total phase shift is 
given by summing the phase shifts in each channel, <5f(/c) = S + (k) + 5~(k). 

Our calculation also requires us to find the bound state energies, which we do by solv- 
ing Eq. (5) using ordinary "shooting" methods. However, once we know the phase shifts, 
Levinson's theorem tells us how many bound states there will be. It works exactly the same 
way as in the bosonic case [4,5]: In the odd parity channel, the number of bound states n_ 
is given by 

(T(0) = 7m_ (10) 
while in the even parity channel the number of bound states n + is given by 

5 + (0) = ttK - \). (11) 

One can derive this result either by the same Jost function methods used in the boson 
case, or by observing that at small k, the nonrelativistic approximation becomes valid so the 
bosonic results carry over directly. For a particular potential there may exist a k = state in 
either of the two channels whose Dirac wavefunction goes to a constant spinor as x — > ±oo. 
(Generically, for k = the components of the Dirac wavefunction go to straight lines as 
x — > ±oo, but not lines with zero slope.) Just as in the bosonic case, such threshold states 
(which we term "half bound states") should be counted with a factor of \ in Levinson's 
theorem. 

Given the phase shifts and bound state energies, we can calculate the one-loop fermionic 
correction to the energy. We work in a simple renormalization scheme in which we add a 
counterterm proportional to 2 — 1, and perform no further renormalizations. The coun- 
terterm is fixed by requiring that the tadpole graph cancel. We must sum the zero-point 
energies — ^uj of the small oscillations, subtract the same sum for the free case, and add the 
contribution of the counterterm. As in Ref. [2-4], we use the density of states 

»«=»<*> + ^ (12) 
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to write 



1 f°° dk d5 F (k) A rr , . 

AH = --^-l _ w _±i + Aftl (13) 

where uj = \fk 2 + m 2 and the sum over j counts bound states with appropriate factors of \ 
as discussed above. Next, we use Eq. (10) and Eq. (11) to rewrite this expression as 

a r T 1^, x f 00 dk , ,dS F (k) . TT .„ „. 

The (divergent) contribution from the tadpole graph is given by the leading Born approx- 
imation to 8 F . Since we have chosen the counterterm to exactly cancel the tadpole graph, 
we can rewrite the counterterm contribution in terms of the Born approximation, giving 

i _ roo Ah A , . 

A/f = --5>i-"0- 1 ^(uj-m)-(5 F (k)-5?(k)) (15) 

where the subtraction S F \k) consists of the first Born approximation to 5 F (k) plus the part 
of the second Born approximation related by <fi — > —(f) symmetry, 

5 {1 \k) = ~J (y F {xf - m 2 ) dx (16) 

which can also be obtained numerically by iterating Eq. (9). As expected, it is indeed 
proportional to (j) 2 — 1. 



III. NUMERICAL RESULTS 

We study a sequence of background fields labeled by a parameter xq that continu- 
ously interpolates from the trivial vacuum <f>(x) = 1 at x = to a widely separated 
soliton/antisoliton pair at x — > oo, 

,/ v ■■ rn(x — x ) . m(x + xn) .„„. 
(j)(x,x ) = 1 + tanh^ ^ - tanh — i ^. (17) 

Fig. 1 shows the fermionic bound state energies as a function of xo- In the limit x — > 
oo two bound states approach energy m^. These are simply the odd and even parity 

combinations of the single soliton bound state at m^. The third (positive parity) bound 
state approaches uj = where the single soliton also has a bound state, but there is only one 
such state (with oo > 0). This is consistent with the result obtained analytically in Ref. [4]. 
As emphasized in Ref. [4], for a single soliton we must thus count the zero mode with a 
factor of |. 

For any finite x , Levinson's theorem holds without subtleties; there are no states that 
require factors of \. Thus for a large but finite x , we find 

* + (0) = I 



5 



1 

0.9 
0.8 
0.7 
0.6 

w. 0.5 

m 

0.4 
0.3 
0.2 
0.1 




^^oo' 

-0 + 

+ 


1 


1 1 


1 1 1 


o 








-\- 














++ +++ ++++ 






















o 








o 






- 


o 





























o 








■ 

I 1 1 









0.5 



1.5 



2.5 

mxo 



3.5 



4.5 



FIG. 1. Fermion bound state energies as a function of mxo in units of m. Even parity states 
are denoted with <3> and odd parity states with +. 



<T(0)=7r (18) 

consistent with having two positive parity and one negative parity bound states (see Fig. 1). 
In the limit x — > oo, an even parity "half-bound" threshold state enters the spectrum at 
uj = m just as in the bosonic case [2]. Also, in this limit, the lowest (positive parity) mode 
approaches u — 0, and is only counted described above. Finally, a negative parity 

mode enters the spectrum from below, also to be weighted by |. 

Thus the net change is to add half of a negative parity state, which via Levinson's 
theorem requires the phase shift S~(0) to jump from n to ^ as x — > oo. This jump occurs 
by the same continuous but nonuniform process as in all cases where a new state gets bound, 
which is illustrated in Fig. 2. Just as in the bosonic case, in the limit of infinite separation 
the potential we have chosen becomes reflectionless, which requires 5 + (k) = 5~{k). 

To contrast the behavior of the zero modes, Fig. 3 shows uj 2 for the bound states of the 
bosonic small oscillations [2] . Because we have chosen the bosonic potential to be consistent 
with super symmetry, the bosonic and fermionic spectra are related. Again as x — > oo the 
bound state energies approach those of the single soliton, and the wavefunctions are formed 
from the odd and even combinations of the wavefunctions for the single soliton. However, 
we see that in the boson case both the mode at m^- and the mode at uj = are doubled, 
so there is no factor of | in counting the bosonic zero modes. For x above approximately 
0.75/m, the system is unstable with respect to small oscillations in one direction in field 
space. This property is reflected in the lowest bound state having uj 2 < [2] . Just as in the 
fermionic case, in the limit xq — > oo, a "half-bound" state enters the spectrum at to — m, and 
the phase shifts in the two channels become equal as the potential becomes reflectionless. 

Fig. 4 shows the values of AH we obtain from Eq. (15) as a function of Xq. In the 
x — > oo limit, we can also do the calculation analytically [4]. We find a phase shift 
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FIG. 3. Boson bound state squared energies as a function of mxo in units of m. Symmetric 
states are denoted with <$> and antisymmetric states with +. 
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FIG. 4. One-loop fermionic correction to the energy as a function of x$ in units of m. 



5 F (k) =4tarT 1 — + 2tan" 1 — (19) 



with Born approximation 



%\k) = f (20) 
and contributions from two bound states at lu = and one at u — m. We thus find 

AH = 2m(- \=] (21) 



7T 



which agrees with the numerical calculation. As shown in Ref. [2], the bosonic contribution 
to the energy is 2m (jt/j — <§^J in this limit, giving a total of — ^ for the soliton/antisoliton 
pair in the supersymmetric model, in agreement with Ref. [4]. 



DISCUSSION 

We have shown how to compute one-loop fermionic contributions to the energies of sym- 
metric background field configurations. Taking these methods together with those of Ref. [2] , 
one can calculate the energy of a localized field configuration in any model of scalars and 
spinors in 1+1 dimensions. This procedure implements standard renormalization conditions 
without ambiguity and lends itself to robust and efficient numerical calculation. A natu- 
ral extension of this work is to treat translationally invariant classical field configurations. 
Long ago, Christ and Lee showed how to quantize the small fluctuations about momentum 
eigenstates formed by taking superpositions of solutions to the classical equations of mo- 
tion at different positions [6]. New problems arise when studying the effects of translation 
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invariance upon fluctuations about time-independent classical field configurations that are 
not solutions to the classical equations of motion. Work is underway to extend the above 
techniques to the nonlocal potentials encountered in this case [7]. 
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APPENDIX: PROPERTIES OF A SINGLE SOLITON 

As a check on the calculations presented in the body of the paper, we present the prop- 
erties (phase shifts, eigenenergies, etc.) for a single soliton directly. In this case, we have 
a potential that satisfies Vf(x) = —Vf(—x) with Vf(x) — > ±m as x — > ±00. We will use 
the bosonic potentials that we obtain by squaring the Dirac equation to obtain a result that 
agrees with the previous method in the limit of a widely separated soliton/antisoliton pair. 
We find 




Mx) = k 2 M%) (22) 



where 



and 



V B (x) = V F (x) 2 - V' F (x) - m 2 (23) 



V B {x) = V(x) 2 + V F (x) - m 2 (24) 

are bosonic potentials. 

Since Vf(x) is antisymmetric, both Vb(x) and Vb(x) are symmetric. We decompose their 
solutions (as bosonic potentials) into symmetric and antisymmetric channels. For x — > ±00, 
we have 

r)l (x) = cos(kx ± 5 s (k)) 
r) A (x) = sin(/cx ± S A (k)) 
Vk( x ) = i cos(kx ± 5 s (k)) 

r}£( x ) = -i sin(kx ± 5 A (k)) (25) 

where the arbitrary factors of ±i are inserted for convenience later. By the (first-order) 
Dirac equation, for all x, 

uj k f k (x) =i(j^ + Vf(x)^ r] A (x) 

^kvt(x) =i(^~ Vf(x)\ v!(x) (26) 
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and 



u k fj A (x) = i(j^ + Vrix)^ v!(x) 



^ S k (x)=i[-^-V F (x))v£(x) (27) 



so that the solutions to the Dirac equation are 
and 



=('|). (29) 

The phase relation between the upper and lower components of these wavefunctions must 
be different as x — ■> ±00 since the mass term has opposite signs in these two limits. 
Putting this all together gives, as x — > ±00, 

cos(A;:r ± 5 s '(k)) = — ( - V F (x)) sin(kx ± 5 A (k)) 
u k \dx J 

= — (k cos(kx ± 5 A (k)) =f m sin(fcc ± 5 A (k)j) 



=F sin(A;a; ± 5 A (k) T 0) = cos(kx ± 5 A {k) T d ±^-) (30) 



and thus 



and similarly 



5 s (k) = 5 A (k)+ tan' 1 ^ (31) 



<^(fc) = 6 s (k)+tajT 1 ^. (32) 

The fermion phase shift in each channel is given by the average of the bosonic phase shifts 
of the two components 

S+(k) = -(S S (k) + S A (k)) = 5 s (k) - - tan" 1 ^ 
2 2 Jx, 



5-{k) = \{8 A {k) + ~5 s (k)) = 5 A (k) - \ tan- 1 ^ 



(33) 



so that 



8 F (k) = 5 + (k) + 5~(k) = 5 s \k) + 5 A (k) - tan- 1 ^. (34) 



k 

We note that in this case there is always a zero-energy bound state given by 
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(35) 



Because it is at zero energy, this state is shared with the antiparticle spectrum and thus only 
contributes with a weight of \ to Levinson's theorem and the sum over small oscillations, 
as we found by following a soliton/antisoliton pair as a function of separation. The extra 
tan -1 y (which goes to | as k — > 0) is thus necessary to connect the fermionic phase shift 
to bosonic phase shifts, which have no such factors of \ associated with their zero modes. 
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